Rashba spin-orbit interaction in graphene and zigzag nano-ribbons 
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We investigate the effects of Rashba spin-orbit (RSO) interactions on the electronic band-structure 
and corresponding wavefunctions of graphene. By exactly solving a tight-binding model Hamiltonian 
we obtain the expected splitting of the bands -due to the SU(2) spin symmetry breaking- that is 
accompanied by the appearance of additional Dirac points. These points are originated by valence- 
conduction band crossings. By introducing a convenient gauge transformation we study a model for 
zigzag nanoribbons with RSO interactions. We show that RSO interactions lift the quasi-degeneracy 
of the edge band while introducing a state-dependent spin separation in real space. Calculation of the 
average magnetization perpendicular to the ribbon plane suggest that RSO could be used to produce 
spin-polarized currents. Comparisons with the intrinsic spin-orbit (I-SO) interaction proposed to 
exist in graphene are also presented. 

PACS numbers: 73.20.At,85.75.-d, 73.63.Bd, 81.05. Uw, 73.43.f 



I. INTRODUCTION 

Understanding the mechanism for the generation and 
manipulation of spin polarized currents is one of the 
greatest challenges for the development of spin-based de- 
vices. Much of the advance in the field in latest yearsi 
has been achieved by studying semiconductor materi- 
als which make up the bulk of current electronic cir- 
cuitry. Among the mechanisms proposed to induce spin- 
polarized currents, the spin Hall effect (SHE) appears 
as the most efficient one. The SHE refers to the phe- 
nomenon in which a spin-polarized current is created 
when an external bias voltage is applied to the system. 
The effect is based on a coupling between spin and mo- 
mentum degrees of freedom, and usually the existence 
of some kind of spin-orbit (SO) interaction in the par- 
ticular system under study is invoked. For bulk semi- 
conductor materials for example, the SO interaction, has 
been proposed to lead to two different manifestations of 
SHE: a) the intrinsic SHE££, in which the material inher- 
its a strong SO interaction from its atomic constituents 
or due to its crystalline symmetries (lack of inversion 
symmetry) and b) the extrinsic SHE^£, in which spin- 
polarized currents appear as a consequence of electron 
scattering by a SO dependent scattering potential. In 
this last situation the scattering potential may be caused 
for instance, by magnetic impurities that couple via a 
spin-orbit term to the conducting electrons or by defects 
that produce spin-dependent scattering. Among these 
scenarios one possibility is when interfaces or surfaces 
are considered. In this case, the existence of the inter- 
face / surface introduces inversion symmetry breaking and 
thus, materials that do not fall into the categories cited 
above can also exhibit SHE. The effective SO interac- 
tion generated in this situation is known as the Rashba 
spin-orbit (RSO) interaction responsible for the Rashba 
effect^. Among the many materials in which RSO inter- 
actions could be exploited to obtain spin-polarized cur- 



rents, graphene presents a unique and intriguing case. 
The material, first isolated as a single layer of graphite 
in 2004 7 . gives access to a crystalline surface with lin- 
ear dispersion around two independent points in its Bril- 
louin zone, the Dirac points. The special dispersion plus 
its crystal structure (two-atom base triangular lattice) 
makes possible to calculate low-energy properties using 
Dirac-type models as the vast literature in recent years 
shows^. Furthermore, the two sublattices of the honey- 
comb structure makes appropriate the use of two- valued 
wave-functions, or spinors, for calculations of various 
properties. The rclativistic description for low energy 
properties has also been used to argue for the relevance 
of additional terms in the standard Dirac Hamiltonian, 
including an intrinsic spin-orbit (I-SO) interaction rep- 
resented by a second-neighbor spin-dependent hopping 
term that respects all the symmetries of the graphene 
planed. One of the consequences of the I-SO interaction 
is to make possible the existence of spin-polarized edge 
states in a new phase of matter, the quantum Spin Hall 
(QSH) phased. In previous works we have studied the 
physics introduced by this interaction in narrow graphene 
ribbons with armchair and zigzag edge terminations and 
in the presence of electron-electron interaction a 11 ' 12 . We 
have shown that the I-SO interaction does not change 
the metallic behavior of armchair nanoribbons in con- 
trast with the predicted result for graphene sheets. More- 
over, the interaction produces spin-filtered states local- 
ized along the edges of the ribbon, (independent of edge 
termination) and, as a consequence, the current induced 
by an applied low external voltage is spin polarized. 
These results are in good agreement with several nu- 
merical and analytic studiesi&Hi&I2iilI&I2i22>2I>22 that 
point to various magnetic instabilities that graphene rib- 
bons may sustain, leading to some kind of magnetic or- 
der along the edges. Unfortunately, all numerical esti- 
mates for the strength of the I-SO interaction remain, 
although still controversial, rather small, in the range of 
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function in momentum space are given by: 



FIG. 1: Left panel: graphene lattice. Lengths along the x and 
y directions are measured in units of a and b = ay/3/2 respec- 
tively. Right panel: First Brillouin zone. K = (47r/3a, 0) and 
K' = (2n/3a,n/b) are Dirac points. 



0.05 - 0.00117716^(600 - 13mK)^. 

However, as a pure two-dimensional material, a 
graphene flake on a substrate lacks inversion symmetry 
and it is natural to expect that a RSO interaction may 
introduce important changes to the material properties. 
The RSO coupling \r is controlled by the applied bias 
and although predicted to be in the range of A_r ~ lmeV 
(Refs i 24 ' 25 i 26 ) recent experiments have shown that it can 
reach values up to Xr ~ 200me V for graphene deposited 
on a Ni substrate 2 ^. Furthermore, experiments also have 
shown short spin relaxation times that suggest an impor- 
tant effect of spin-orbit interactions in graphen e 28 ' 29 i 30 ' 31 . 

These new developments highlight the need for a better 
understanding of the role played by the RSO interaction 
on various properties of graphene and graphene ribbons. 
In this work we address the questions raised by the pres- 
ence of RSO in graphene sheets and zigzag ribbons. The 
RSO interaction strongly affects the dispersion relation 
near the two-independent Dirac points as well as the na- 
ture of the corresponding wave-functions as we will show 
below. As a SU(2) breaking symmetry interaction, it fa- 
vors a spatial spin ordering but in contrast to the I-SO 
interaction introduced above, the spin order is not orig- 
inated on each individual state having the same spatial 
spin distribution but it emerges from averaging over sev- 
eral states. 



II. GRAPHENE SHEET: MODEL 

To describe a graphene sheet (an infinite mono-layer 
of carbon atoms arranged in a honeycomb structure) , we 
introduce, as usual, two sublattices A and B with their 
respective atoms connected by vectors: 

5 1 =a(0,f/V3) 

5 2 = a(l/2,-l/2 v / 3) 

« 3 = a(-l/2,-l/2V5) (1) 

where a = 2.4 A is the lattice constant [see Fig. (JTJ)]. 

In the absence of spin-orbit interactions spin-up and 
spin-down electrons are degenerate. The SU(2) sym- 
metric Hamiltonian and four-component spinor wave- 
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with ip(k x , fey) = /j(e lfc « 26 / 3 + 2 cos l^ e -ik v b/3y In these 
expressions (p is defined as <p(k x ,k y ) = ip(k x ,—k y ). No- 
tice that for real values of k y , Cp = ip* . The eigenvalues 
of Eq.([J]) are E = ±e = ±y/ipip and the corresponding 
eigenvectors arc defined in terms of the angle ao as: 
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with (p = |cp|e m ° and N the normalization factor. For 
neutral graphene, ^ + (^-) represents solutions with 
E > (E < 0) and refers to electron (hole) conduc- 
tion (valence) bands. In this language, the particle-hole 
symmetry implies that for each electron state with en- 
ergy E — e and eigenstate characterized by ao, there is a 
hole state with E = —e and eigenstate given by ao + ir. 



A. Rashba spin-orbit interaction 

Depositing graphene on substrates and/or applying 
external fields makes possible to introduce a control- 
lable RSO interaction. In the following we take the 
effective electric field £ perpendicular to the graphene 
plan o 10 ' 32 i 33 . The Rashba Hamiltonian is then given by: 



H R = ic l(uv ■ cOcj + h.c 

<ij> 



(5) 



where a represents the Pauli spin operator for the spin 
degree of freedom. Here Uij is given by: 



^ = 2^fe^*^ = -T***«- (6) 



where £ is the applied electric field in the direction per- 
pendicular to the graphene sheet, d = a/y/3 is the dis- 
tance between the two adjacent sites and Sij = 
dj — Si is a vector on the graphene plane. 

With these definitions the RSO interaction takes the 
form: 



H R = ic\ a R a(i c B ii + h.c 



(7) 
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where 

+e - 8 fc x dv^/2- i fc y d/2(_ CT ?,y 2 + f 2 ) (8) 

where (£,?7) stand for spin up and down. 

The RSO interaction couples spin-up and spin- 
down states, breaking the corresponding SU(2) symme- 
try, leading to the Hamiltonian (written in the four- 
component spinor $ basis): 
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where 

<p = te l2fc * b / 3 (l + 2e^ fc « b cos(fc :!; a/2)) 

^+ = Aije 12 ^ 6 /^! + 2cos(fc x a/2 + 27r/3)e" lfc " b ) 

V- = A fl e 42fc * b/3 (1 + 2 cos(k x a/2 - 27r/3) e - 4fc « b $10) 

The eigenvalue equation is given by: 



FIG. 2: Energy bands as function of k x of an infinite graphene 
plane without (small dotted lines) and with (filled circles and 
empty diamonds) RSO. The strength of the RSO interaction 
is \r = 0.2t. The plotted bands correspond to k y = n (lower 
set) and k y = tt/2 (upper set). The inset shows the new lower 
energy band with k y — tt and the splited Dirac point. A finite 
gap separates all other conduction and valence bands. 
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in which we have defined 



(11) 



point at K = (27r/3;7r/6) shown in Fig. @, the position 
of one new point is at (k' x , n/b) with k' x given by: 



El 



E 2 



$ = itp+tp- i<p<p- = E + E~e iu (12) 
and £o;± = |V0;±|- The angle v is defined by 

taniv = S$/5R$. (13) 

Note that v is independent of the RSO coupling. 

Equation (fT2"]) shows explicitly that E — ► — E, repre- 
senting the particle-hole symmetry, is preserved by the 
RSO interaction. Using this property, in the rest of the 
paper we will focus on conduction bands only. 

In Fig. ([2]) and Fig. ([3]) we plot few bands for the 
infinite graphene plane in the presence of the RSO in- 
teraction. Two new features appear: a) due to the 
breaking of the SU(2) symmetry, each of the degener- 
ate bands (in the absence of spin orbit interaction) splits 
into two. For a given (k x , k y ) the energies of the newly 
separated upper (E u ) and lower bands (E[) are related 



by Ei 
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ei). b) The inset of 



Fig. ([2]) shows the change in the band with k y = tt which 
originally touches the corresponding valence band at a 
Dirac point located at K = (27r/3a, ir/b). 

A remarkable consequence of the RSO interactions is 
the splitting of the original Dirac point caused by cross- 
ings of conduction and valence bands. The location of 
the new points in reciprocal space respects the underly- 
ing honeycomb symmetry and depends on the strength 
of the interaction A# as shown in Fig. ([J|. For the Dirac 
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In the linear (Dirac) approximation of the Hamiltonian 
and for small values of the interaction strength A#, the 
splitting of the Dirac points is missed, and the low-energy 
effective Hamiltonian describes graphene with RSO inter- 
actions as a zero-gap semiconductor as reported in pre- 
vious works^. 

It is important to remark that the RSO interaction 
does not open a gap in the spectrum at the Dirac point, 
in contrast to the I-SO interaction mentioned above. 

To solve for the eigenstates of the Hamiltonian ([9]) we 
notice first that in the limit A# — > the spinor introduced 
in Eq.© takes the form * = a iv l 2 ^ x + e~ ll/ ^ 2 ip l where 
ipH are the two degenerate spinors defined in the absence 
of the RSO interaction. The four components of & satisfy 
the following relations: 
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FIG. 3: Energy bands as function of k y of an infinite 
graphene plane without (dotted lines) and with (full lines). 
The strength of the RSO interaction is = 0.2t. The plot- 
ted bands correspond to k x = 0.67T. As described in the text 



there are four degenerate states at 
for each k x value. 



±fci and k y = ±&2 



- 




FIG. 4: The positions of new zero-energy (Dirac) points in 
momentum space for different values of the RSO interaction 
strength \r. The distribution of the new Dirac points around 
the original ones has the 27r/3 rotational symmetry of the 
graphene lattice. 



which lead to the cigenstates: 
/ 





(16) 







a! 





e — • — e— • — e- 



=W/2 



y=-W/2 



FIG. 5: Hard-wall boundary conditions for ZGR imposes 
ua = on the lower edge and ub = just before the upper 
edge (green dotted line). The deformed lattice shown on the 
right side, corresponds to a gauge transformation (see text) 
and it is equivalent of labeling both A and B sites in each 
zigzag line with the same y-coordinate. The horizontal oval 
shows our choice of the unite cell and the vertical oval shows 
the choice of the unit cell which corresponds to Eq. (|20[1 . 



Here a is defined by: 



EE- E 
E E+ E 
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In the limit \r — > 0, y — * 1 and a — ► cto- The 
expressions above (Eq. (|16jl ) correspond to a state with 
energy E in the split upper band and momentum (k x , k y ). 
The sate with the same momentum components in the 
split lower band is obtained by the replacement v — * ^+7r. 
The remaining particle-hole symmetric states (in the split 
valence bands) are obtained by taking a — > a + n. 

Wave function (|16p has the property of ua] — Ubi 
and ubi = waj.- This reflects that, in the presence of 
the SU(2) symmetry breaking RSO interaction the prob- 
ability of finding an electron in the spin up state i.e 
|?Mt| 2 + I m st| 2 is equal to the probability of finding it 
in the spin down state i.e |maj,| 2 + |w_bJ 2 . This is a direct 
consequence of the fact that the RSO interaction does not 
break time reversal symmetry. However, as we will sec 
below, this dos not exclude the possibility of separating 
the spin-up and spin-down electron states and localizing 
them at different positions in the sample while preserving 
a zero net magnetization. 



III. ZIGZAG GRAPHENE NANORIBBONS 

A. Zigzag nanoribbons without Rashba spin-orbit 
interaction 

To study the interplay between confinement and the 
RSO interaction we analyze the case of zigzag graphene 
(ZGR) nanoribbons, defined according to Fig. ([5]). 

As standard practice the hard-wall boundary condi- 
tions are imposed by setting Ua(jj — 0) = on the lower 
border and ub(W — b/3) = on the line located at a dis- 
tance 6/3 below the upper borde r 34 ! 35 ^^^. After close 
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inspection, however, one realizes that it is more conve- 
nient to label all the atoms alone each zigzag line with 
the same y-coordinate. This is equivalent to work with 
the deformed lattice shown in the right panel of Fig. ([5]) 
and it amounts to perform a global gauge transformation 
CB(ky) — > CB{ky)e 2lkyb ^ on the original Hamiltonian. It 
also represents a different choice of unit cell as described 
in Fig. ([5]). For symmetry reasons we also set the y-axis 
to be in the center of the ribbon. With this choice, the 
hopping term <p reduces to: 



ip(k x , k y ) = i(e 



2 cos ) 

2 ' 



and the boundary conditions are: 

u A (y = -W/2) = 0, u B (y = W/2) = 0. 



(18) 



(19) 



Notice that in most of the literature on graphenc rib- 
bons the usual convention for tp and the boundary con- 
ditions are: 



u A (y = -W/2) = 0, u B (y = W/2 

v (k X7 k v )=t(l + 2cos^e-* k «») 



1)=0 



(20) 



which correspond to choosing a unit cell along the vertical 
link in the right-side panel of Fig. ([5]). 

The wave-function of the ZGR can be found in a 
straightforward manner as follows. Since k x is a good 
quantum number, the wave-function for a given k x must 
be a superposition of degenerate states with different k y 
values. In the absence of SO there are only two degener- 



ate spinors for each k x namely at k y = k and k y 



-k. 



Therefore the wave-function is the superposition of these 
two spinors: ^ = a^(k x ,k) + b^(k x ,— k). After apply- 
ing the boundary conditions as given in Eq. (|19p . we find 
b = —a such that 



Ne 



sin(ao/2 + ky — mr/2) 
sin(— Qo/2 + ky — mr/2) 



where k satisfies 



kW 



(21) 



(22) 



Figure © shows the conduction bands of a ribbon with 
W = 4b. Zigzag ribbons present two remarkable fea- 
tures as compared to graphene sheets: the momentum 
across the ribbon k y can take complex values between two 
Dirac point a 34 i 35 ' 36 i 37 ' 38 , producing edge states and their 
band-structure depends on the width W or the number 
of zigzags chains N = W/b — 1. It can be show n 39 i 40 ' 41 i 42 
that in zigzag ribbons with odd number of chains N, the 
so-called 'zigzag/zigzag' configuration, conduction and 
valence edge bands cross a = 7T. In contrast, ribbons 
with even number of chains N, in the 'zigzag/anti-zigzag' 
configuration, edge bands do not cross albeit remain de- 
generate at k x a = it. 
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FIG. 6: Energy bands of a zigzag ribbon with W = 46 in 
the absence of the RSO interaction. Each band is doubly 
degenerate due to the SU(2) spin symmetry. The edge band 
(diamonds) corresponds to an imaginary value for the label k 
in Eq. (p2)l . 



B. Zigzag nanoribbons with Rashba spin-orbit 
interaction 



As seen in the previous section, ZGRs present the pe- 
culiar feature of edge states which remain highly quasi- 
degenerate at low-energies for wide ribbons. These states 
are expected to be strongly affected by the presence of 
a RSO interaction. Below we proceed to obtain the ex- 
act expressions for the band-structure and corresponding 
eigenstates for the ZGR with RSO interactions. 

It is necessary to remark first that, since the RSO in- 
teraction involves nearest neighbor hopping, the bound- 
ary conditions as imposed in Eq. (|19|) remain unchanged. 
However, for a given value of k x , there are four degener- 
ate states at k y = ±k± and k y = ±&2 in contrast with 
the previous case (with only two degenerate states at 
±fc). This is easily seen in Fig. ([3]). Notice that there are 
certain energies, such that it seems that there are only 
two degenerate states, however the wave-function is re- 
ally the superposition of four spinors with ±fci and ±&2 
taking imaginary or complex values. The general wave 
function is: 



ZGR 



aiPih) + ty(-fci) + CT/>(fc 2 ) + d^(-fca) (23) 



where tp(ki) = ip(k x , ki) and k\ and ki satisfy the condi- 
tion given by the degeneracy: 



E(±ki) = E{±k 2 ). 



(24) 



Imposing the boundary conditions given in Eq. |T 
yields 



sm(/ij + S) _ E} sin(/4 + S) 
Ei sin(/x 1 " + 5) E2 sin(/i2 + 5) 



(25) 
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FIG. 7: Energy bands of a ZGR with W = 46 and X R = 
0.4i. Different regions correspond to values of k\ and fo: I) 
both real, II) one real and the other one imaginary III) both 
imaginary, IV) complex with ki = k%. 



where fif = {vi ± cti T- kiW)/2 and 6 = ±7r/2. These two 
equations define the band structure and the correspond- 
ing wave-function in terms of the width W and the RSO 
coupling A#. The wave function coefficients are given by: 



o = -b = N\l sin(/ij + 5) 



d = -c = N\-±;Bm{nZ + 5) (26) 
E i 




sions it can be shown that ^^ R (y) 



where N is the normalization factor. From these cxprcs- 

*ff*(-I/) and 

Figure J?]) shows the conduction bands for a ribbon 
with W = 46 and Ah = 0.4t. Starting at k x = both 
parameters fci, fc 2 are real (region 7). As fc^ is increased, 
/c2 goes to 7T or zero and in region II it becomes complex 
(with constant real part equal to w) or purely imaginary. 
In region III both k\ and k% take imaginary values with 
a constant real part of 7r or zero. The energy of the lower 
band goes to zero E = at the point defined by: 



a»(-2-) 



4A| 



4t 2 



(27) 



Notice finally that there is also a region (IF) where fci 
and ki are complex conjugate of each other. As it occurs 
with the I-SO interaction, the presence of the RSO in- 
teraction lifts the apparent quasi-degencracy of the edge 
band while preserving the Dirac points 12 . The expression 
for the dispersion of the edge bands of an TV- wide ribbon 
is readily obtained and is given by: 



LU 
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FIG. 8: The upper panel shows the edge bands of a ribbon 
with W = 36 (N = 2) and Xr — 0.2t. Conduction and valence 
edge bands do not cross at the band center. In contrast, the 
lower panel shows a W — 46 (JV = 3) ribbon with crossing 
bands. The full lines are fits using the expression in Eq. (|28[) . 



E « ±t(k x a - k%a 



0„\N 



(28) 



where is defined in Eq. (|2"7| . 

It is interesting to notice that RSO interactions pre- 
serve the power law energy dispersion and edge bands 
crossing/anticrossing feature as shown in Figs. (0). 

With the expressions obtained for the wave-functions, 
it is possible to calculate various quantities. In particular, 
Fig. ([9]) shows the spatial probability distribution for S z , 
the z-component of the spin operator defined as < S z >— 
|u^|| 2 + |u_bt| 2 — l M Ai| 2 — l u si| 2 f° r the lowest energy 
conduction band of a ribbon with W = 46 and \r = 0.4t. 

The figure highlights the fact that the RSO interaction 
produces a clear spin polarization on the edge states of 
the ZGR. The non-homogeneous spin distribution across 
the ribbon is, however, highly dependent on the state 
considered. This is in contrast to the effect produced 
by the I-SO interaction where each state becomes spin- 
polarized with the same spatial spin distribution^. 



7 




-0.6 1 1 1 1 1 1 1 1 1 

-2 -1.5 -1 -0.5 0.5 1 1.5 2 

y/b 



FIG. 9: Expectation value of the z— projection of the spin 
operator, < S 2 > as a function of position across the ribbon. 
The curves are calculated using the lowest energy conduction 
band of a ZGR with W — Ab and \r = 0.4t. Different curves 
correspond to different values of k x of Fig. |(7J 

IV. CONCLUSIONS 

Grapheme ribbons show unique and interesting trans- 
mission properties due to its band-structure and the 
pseudo-spin nature of its wave-function o 39 ' 40 . The rel- 
ativistic nature of the description normally used makes 
it necessary to understand further other relativistic ef- 
fects that could alter their transport properties. In this 



work we have investigated the consequences of one of such 
interactions: the Rashba spin-orbit interaction that is 
expected to be relevant under applied external bias volt- 
ages. We have shown that in graphene sheets, the RSO 
removes the SU(2) spin degeneracy as expected while it 
does not open a gap in the spectrum. It does, however, 
introduce additional Dirac points in the Fermi surface at 
low energies due to crossings between valence and con- 
duction bands. 

Because of its peculiar edge band, zigzag graphene rib- 
bons are potential candidates for spintronic applications. 
The edge bands are expected to be magnetically unsta- 
ble and as such to be strongly affected by SO interac- 
tions. We have shown that the RSO in particular pro- 
duces states that have spin polarization and are strongly 
localized along the edges. These states present opposite 
polarization at opposite edges and the spatial spin distri- 
bution is strongly dependent on the state under consid- 
eration. Without external fields the net spin polarization 
of the ribbon remains null as a natural consequence of the 
conservation of time reversal symmetry under the RSO 
interaction. However, these results suggest the possibility 
to obtain spin polarized currents if the states selected by 
an applied external voltage sustain an average non-zero 
spin polarization. 
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